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ABSTRACT 

We have searched for [OIII] 5007 emission in high resolution spectroscopic data from Flames/Giraffe 
VLT observations of 174 massive globular clusters (GCs) in NGC 4472. No planetary nebulae (PNe) 
are observed in these clusters, constraining the number of PNe per bolometric luminosity, a < 0.8 x 
1O~ 7 PN/L0. This is significantly lower than the rate predicted from stellar evolution, if all stars 
produce PNe. Comparing our results to populations of PNe in galaxies, we find most galaxies have 
a higher a than these GCs (more PNe per bolometric luminosity - though some massive early-type 
galaxies do have similarly low a). The low a required in these GCs suggests that the number of PNe 
per bolometric luminosity does not increase strongly with decreasing mass or metallicity of the stellar 
population. We find no evidence for correlations between the presence of known GC PNe and either 
the presence of low mass X-ray binaries (LMXBs) or the stellar interaction rates in the GCs. This, 
and the low a observed, suggests that the formation of PNe may not be enhanced in tight binary 
systems. These data do identify one [OIII] emission feature, this is the (previously published) broad 
[OIII] emission from the cluster RZ 2109. This emission is thought to originate from the LMXB in 
this cluster, which is accreting at super-Eddington rates. The absence of any similar [OIII] emission 
from the other clusters favors the hypothesis that this source is a black hole LMXB, rather than a 
neutron star LMXB with significant geometric beaming of its X-ray emission. 

Subject headings: globular clusters: general - stars: AGB and post-AGB - galaxies: stellar content - 
stars: evolution 



1. INTRODUCTION 

Emission lines in the spectra of extragalactic globular 
clusters (GCs) provide a useful method for identifying 
exotic objects in their integrated light. In particular, 
the presence of the [OIII] 5007 emission line in a GC is 
indicative of the presence of objects such as planetary 
nebulae (PNe), supernova remnants or X-ray binaries. 
Despite our interest in identifying and studying such ob- 
jects, and the frequency with which they occur, few sys- 
tematic searches for emission lines in the spectra of GCs 
currently exist. 

One of the main sources of [OIII] emission in a stellar 
population are PNe. These nebulae are produced near 
the end of an asymptotic giant branch (AGB) star's life 
when gas, driven outward by a wind or superwind, is 
ionized by the hot inner core. Our understanding of the 
formation of these systems remains uncertain. In partic- 
ular, challenges remain i n explaining the axisymmetri c 
shape of some PNe (e.g. [Sokerl[200a iDe Marcol[2009T ). 
the universality of the brightest PNe o bserved across a 
wide range of galaxy morphologies (e.g. [Qiardull o et al.1 
120021: lCiardulloll2010D and the number density of PNe in - 
different galaxies (iHui et alJ ll993t iCiardullo efall 120051 : 
iBuzzoni etafl 120061 ). In recent years, theories for the 
formation of PNe have suggested that stellar interac- 
tions in binary systems may play an important role (e.g. 
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IMoe fc De Marcol [20061 : ISokerl 120061 : IDe Marcol l2009h . 
The population of PNe formed in GCs provides an in- 
teresting addition to studies of PNe. GCs provide a rel- 
atively well constrained stellar population in which to 
study PNe in old, metal poor, environments. They also 
provide a unique dynamical environment in which to con- 
sider the effects of binaries and stellar interactions on the 
formation of PNe. 

In the Galactic GC system, searches for P Ne have iden- 
tified only four such systems in 133 clusters (pacobv et al.1 
19971). Based on the total mass of these clusters, 



Jacobv et al.l fl997) proposed that PNe may be under 



abundant in globular clusters, compared with galactic 
fields (although at a significance of only around 3<r). 
In extragalactic GC systems, PNe can be identified in 
their unresolved spectra via the strong [OIII] emission 
from these nebulae. However, there have been few stud- 
ies dedicated to identifying emission lines from GCs and 
currently only a few extragalactic GC PNe have been 
recorded. These were identified via the 5007A emis- 
sion line from GC s in NGC 3379 (iBergond et all [20061: 
iPierce et af1l2006D . the Fornax dwarf (|Larsenl 120081 ) and 
NGC 7457 (which m ay be an intermediate age cluster, 
iChomiuk et all [20081/ In this paper, we examine high 
resolution spectra of GCs around NGC 4472, with the 
aim of investigating the formation of PNe in GCs. These 
clusters have a combined luminosity of ^2x 1O 8 L0 j v, 8x 
more than the Galactic GC system. 
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Another potential source of [OIII] emission from GCs 
are nebulae produced around bright X-ray sources. X- 
ray and [OIII] emission has previously b een observed 
from a GC in NGC 4472 flZepfet all [2001 and a GC in 
NGC 1399 (jlrwin et all 120101) . The [QUI] emission ob- 
served from the NGC 447 2 GC is very bro ad, implying 
velocities of -2000 km/s (|Zepf et al.ll2008[) . This is too 
high to be from a PN in the cluster. Instead it is likely 
that this emission is associated with the bright black hole 
low mass X-ray binary (LMXB) in the cluster. This X- 
ray binary is thought to have super-Eddington accretion 
rates and is likely ionizing a wind driven fro m the sys- 
tem (|Maccarone et all 120071: iZepf et alJfeOOTl) . The de- 
tection of similar (broad) [OIII] emission in GCs provides 
a method for identifying more of these interesting sys- 
tems and investigating the properties of X-ray binaries 
with very high accretion rates. 

2. NGC 4472 DATA/ REDUCTIONS 

This study is based on archival spectra of GCs in the 
galaxy NGC 4472. The survey targeted known clusters 
from the catalog of (|Rhode fc Z epf 2001) in four differ- 
ent fields around NGC 4472. Details of each observed 
field are given in table [1] These data were obtained 
on various nights between May to July 2003 and Jan- 
uary to March 2006, under the programs 072.B-0384(A) 
and 075.B-0516(A) (PI: Bergond). Here, we review these 
data, but f urther details o f this a nd similar survey s can 
be found in lBergond eTall (|2TM ) : [ZepfeTall ([20071 ) and 
Zepf et al. (in prep). 

These data were obtained on the very large tele- 
scope (VLT) using the F LAMES instrument i n its GI- 
RAFFE/MEDUSA mode (|Pasauini et al.ll2002D . This is 
a fiber fed, multi-object spectrograph, consisting of 130 
fibers which can be positioned over a 25' field. The spec- 
tra cover the range 5010< A <5810A with 0.2A pixels 
and a resolution of 0.9A. 

We obtained the raw science data from the ESO 
archival and reduced it using the GASGANO software^- 
Master dark and bias frames were produced for each ob- 
servation using the routines gimasterdark and gimaster- 
bias. In the absence of the original reduced calibration 
files, we used the standard calibration data sets pro- 
vided by the GIRAFFE pipeline (as recommended in the 
pipeline documentation). These include the flatfielding/ 
transmission correction, fiber location/ width determina- 
tion and dispersion solution. The data were then reduced 
using the giscience routine in GASGANO. This routine 
performs the bias, dark and flatfield corrections. It then 
applies the dispersion solution and extracts each spec- 
trum from the frame. The spectra were then associated 
with the targets using the fiber association data table 
(included in the ESO archive data). 

The sky was measured independently for each expo- 
sure via dedicated sky fibers, spread across the field of 
view. The sky level was not found to correlate strongly 
with either location in the field of view or distance to the 
center of NGC 4472. The lack of a correlation with galac- 
tocentric radius is unsurprising. This is because GCs in 
the central regions of NGC 4472 were not targeted by 
this survey, therefore the galaxy surface brightness is rel- 

I http://archive.eso.org 

II http:/ /www. eso.org/sci/software/gasgano/ 



TABLE 1 
Targetted fields 



Field 


RA 1 


DEC 1 


Number 


Total exp 


Number 








exp 


time (ks) 


clusters 2 


NE 


187.56317 


8.10417 


4 


10.0 


108 ( 80) 


SE 


187.55866 


7.88373 


6 


15.6 


114 (101) 


SW 


187.33344 


7.88362 


5 


11.5 


112 ( 74) 


NW 


187.35649 


8.14778 


4 


10.4 


117 ( 87) 



1 approximate center of each field. 

2 parentheses indicates the number of clusters for which their 
velocities can be measured from these data (Zepf et al., in prep). 



atively faint and fainter than the sky for most regions. 
For the most central GC, which is — 3' from the center of 
NGC 4472, we expect the galaxy to contribute —25% of 
the GC flux. Given the lack of correlations, the improved 
S/N produced by a simple average of all sky fibers was 
found to produce the best results. The cluster spectra 
were sky subtracted and combined using the IRAF tasks 
IMARITH/SCOMBINE. An approximate flux calibra- 
tion was then applied to each clus ter spectrum based o n 
its known V-band mag nitude from lRhode fc Zepfl (pOOlh . 

This provided 453 spectra of NGC 4472's GCs. How- 
ever, several clusters were observed multiple times (in 
different fields). In total, these data provide spectra for 
358 unique clusters. Where multiple spectra were avail- 
able for a cluster, these were averaged together. 

The GCs in galaxies like NGC 4 472 cover a wide 
range of radial velocities. For example, ICote et aLl (2003) 
show that NGC 4472's GCs have velocities in the range 
250 < v < 1950 km/s, corresponding to wavelength shifts 
of 4 < AA < 32A. It is therefore difficult to identify 
5007A emission from clusters without determining their 
redshift. As part of a separate study, these data have 
been used to derive the radial velocities of 174 of these 
clusters (Zepf et al., in prep.). For the remaining clusters, 
the S/N ratio was too low to obtain reliable velocities. 
The high resolution of these data allowed velocities to be 
estimated with errors <20 km/s. We used these veloci- 
ties to estimate the rest frame wavelengths of these data 
with AA < 0.3A. 

3. 5007A EMISSION FROM NGC 4472'S GCS 

We used these data to identify 5007A emission from 
the 174 GCs with known velocities. This was done by 
measuring the flux at the (redshift corrected) wavelength 
of 5007A using two boxcars with width 2 and 5A. For 
each cluster, we fit the continuum level via a linear fit 
to the spectrum, from 5010-5060A, excluding the region 
4A either side of the line. The 2A width measurement 
is appropriate for measuring emission from PNe, which 
typically have velocities of — 50km/s. Therefore the large 
majority of their flux should be in this narrow region. 
The wider (5A) bin has greater scatter, due to the in- 
creased noise associated with the wider region, but is 
more sensitive to broad emission features from non-PNe 
sources, such as that seen from RZ 2109. 

Figure [1] shows the 5007A luminosities of these clus- 
ters. In this figure, we present line strengths as lumi- 
nosities, rather than equivalent widths. This is because 
the clusters in our sample have different luminosities (and 
hence different continuum levels). The line luminosities 
are relative to the continuum level, such that emission 
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Fig. 1. — NGC 4472 GC luminosities in the [OIIIJ5007 emission line as measured with respect to the continuum, in a 2A (top) and 
5A (bottom) region centered on rest-frame 5007A. Positive luminosities indicate emission features above the continuum, while negative 
luminosities represent absorption features. The dashed line corresponds to the luminosity of a PN 2.5 mag beneath the brightest PN cutoff 
luminosity. The right hand panels show histograms of the measured luminosities. They demonstrate that the measured luminosities are 
consistent with Gaussian noise distributed around zero. One cluster can be identified in these data as having excess emission over these 
regions. This is the (previously identified) cluster RZ 2109. It can be seen that the luminosity of this emission feature increases significantly 
in the wider line region, suggesting a very broad emission feature. 

and absorption lines have positive and negative fluxes, 
respectively. The right panels of figure [T] show the distri- 
bution of 5007A line luminosities. It can be seen that the 
measured line luminosities have an approximately Gaus- 
sian shape centered close to zero luminosity and with 
a standard deviation of 8.5x 10 34 erg/s/A. This is sim- 
ilar to that expected from the noise of the average re- 
duced cluster spectra. These have a ~ 6 x 10 34 erg/s/A, 
and thus a in our 2A bins is v2 times this value, or 
a(2k) = 8.5 x 10 34 erg/s/A. Such a distribution would be 
expected due to noise, if no 5007A features are present 
in the clusters. We therefore assume that the scatter 
observed in figure Q] is due primarily to noise. Stud- 
ies of PNe often quote the population of PNe within 
2.5 mag of the brightest PN cutoff magnitude (M*), 
where (Mf om i-2.5) = -1.98 (corresponding to a luminos- 
ity, L[om] ~ 2.4 x 10 35 erg/s). This limit is plotted in 
figure [T] as the dashed line. It can be seen that PNe 
above this limit would be detected by these data. 

The top panel of figure[T]idcntifies only one cluster with 
emission above this limit. The luminosity of this cluster 
increases significantly in the wider 5A region, suggesting 
that this emission is very broad. This emissio n has previ- 
ously been identified in this cluster (RZ 2109, Zepf et al. 



and was studied in greater detail by iZepf et all 
), who identified [OIII] 5007 and [OIII] 4959 lines 
with widths of ^2000km/s. No clusters with similar, 
broad, emission are identified in these data. In section 
[5l we discus this result in the context of the population 
of ultraluminous LMXBs in GCs. 

No other emission features are observed in these data 
(we confirmed our automated approach by visually ex- 
amining the cluster spectra for emission lines) . This sug- 
gests that no PNe, in the brightest 2.5 mag of the plan- 
etary nebula luminosity function (PNLF), are present in 
these clusters. The implications of this non-detection on 
the GC PN population is discussed in the following sec- 
tion. 

4. PLANETARY NEBULAE IN GLOBULAR CLUSTERS 

4.1. How many PNe are expected? 

The absence of any PNe in these data raises the ques- 
tion of how many are expected. To investigate this, 
we consider some simple theoretical predictions for the 
number of PNe in a stellar population, per bolometric 
solar luminosity, o cHL An estimate for a can be ob- 

111 We use a to represent the total population of PNe and a x to 
represent the population with [OIII] magnitudes in the range M* 
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tained from sim ple stellar evolution to be, a ~ 2 x 
IQ-^VpnPN/L^ (|P^nzin[ fc Buzzonilll986HJacobv et al.l 
119971: iBuzzoni etldl 120061) . This is the number of post 
AGB stars formed (per solar luminosity) x the life- 
time of an observable PN, tpn- While the lifetime of 
PNe may vary in different stellar populations, we can 
estimate an upper limit as the dynamical timescale, 
7"PN,dyn ~ 30,000 yr. This is the time it takes the 
nebula gas to disperse (and hence end the observable 
PN phase). These considerations suggest that, if all 
stars produce PNe with lifetimes close to the dynamical 
timescale, a ~ 6x 10 _7 PN/L Q . Some stellar populations 
may have shorter PNe lifetimes, or may not produce ob- 
servable PNe at all, therefore this value is an estimate 
for the maximum number of PNe expected, per solar lu- 
minosity. 

4.1.1. tpn and age 

The age of a s tellar population can have a signifi- 
cant effect on a. Bu zzoni et al.l (|2006l ) considered the 
effect, on a, of the decreasing mass of PNe central stars, 
with the increasing age of the stellar population. They 
demonstrated that initially tpn (and hence a) should in- 
crease, with increasing age, due to the slower evolution 
(and hence longer lifetimes) of the lower mass central 
stars. However, as the age increases further, it takes the 
central stars longer to reach the temperatures required 
to ionize the PN, hence reducing their observable life- 
time (and reducing a). Eventually, one would expect 
the central star mass to fall beneath a critical limit of 
~0.52M Q . At these low masses, the star will evolve too 
slowly to reach the temperatures required to ionize the 
circumstellar gas, and hence produce an observable neb- 
ula. This limit may be important for very old stellar 
populati ons, such as GCs and th e oldest populations in 
galaxies. iMoe fc De Marcol (|2006() investigated the num- 
ber of PNe expected in a stellar population with an age 
of 11.5 Gyr and a turnoff mass of 0.85Mq -typical of a 
GC. Their conclusion was that no PNe are likely to be 
produced in such a population, based on the evolution 
of isolated stars. However, they note that there is an 
uncertainty in this prediction, due to the uncertainty in 
the central star low mass cutoff (and the relationship be- 
tween main sequence turnoff mass and the resulting PN 
central star mass). 

4.1.2. Binary interactions 

The above estimates for a consider the evolution of an 
isolated star. However, it is likely that binary interac- 
tions play an important role in producing PNe. It has 
been proposed that most (if not all) of the PNe that are 
observed are pr oduced through some form of binary in- 
teractions le.g.lBondl2000tlDe Marco et all 120041 : ISokerl 
I20M iMarco fc Moell2005l : be Marcoll2009ft There are a 
few reasons to suspect that binary interactions are im- 
portant in the formation of PNe. Firstly, the brightest 
PNe (M*) observed in galaxies across all morphologies 
are found to be very similar, regardless of the age of the 
stellar population. This is not expected from the evo- 
lution of a single st ar. For example, the simulations of 
iMarigo et al.l (2004) suggest that M* should decrease by 

to M*-x 



approximately 3 mag in the first 6 Gyr of star formation 
(|Ciardullol2010D . Secondly, most PNe observed are found 
to be non-spherical. Both of these observational con- 
straints are hard to produce from the evolution of a single 
star, but ca n potentially be explained in binary scenar- 
ios (s ee e.g. ICiardullo et~aDl2005t lSokerll2006l : iDe Marcol 
2009) . If binary interactions are required to produce ob- 
servable PNe, then we require that the PN predecessor 
needs to be in a binary system which is tight enough 
that interactions can occur during the AGB phase. Ad- 
ditionally, the binary separation must be large enough 
that interactions do not occur at earlier times, such as 
during the red giant branch ph ase, when the core is not 
hot enough to ionize a nebula. IDe Marco fc Moel (2005) 
propose that, this limits PN forming systems to binaries 
with separa tions (a) o f rough ly 100 < a < 500R Q . In a 
later paper. IDe Marcol ()2009h estimate that 13% of stars 
reside in binaries in which such interactions can occur, 
with a corresponding reduction in a. The binary for- 
mation mechanism therefore predicts significantly lower 
a than would be expected based solely on stellar evolu- 
tion. 

4.1.3. a in galaxies 

Empirical estimates for a are also available for several 
nearby galaxies. However, it should be noted that, ob- 
taining accurate measurements for a is challenging. In 
the Milky Way, estimates for a vary significantly, par- 
tially due to difficulties in measuring the d istances to 
Galactic PNe. For example, iPhillipsI (|2002t) est i mated 
a ~ 6.3 x 10~ 7 in the Milky Way, while ISokerl (pOOfih 
suggested tha t the distances to the PNe considered by 
Phillips (2002) may be, on average, 1.5 x greater, imply- 
ing a ~ 1.9 x 10~ 7 . The best constrained PNe popula- 
tions are those of the SMC and LMC. For these galaxies 
large numbers of PNe (with known distances) can be de- 
tected across a significant fraction of the PNLF. For more 
distant galaxies, beyond our local group, PNe surveys are 
often constrained to the brightest tip of the PNLF (with 
detection limits ranging from 0.5-2.5 mag fainter than 
M*). Their total population of PNe is then extrapo - 
lated, based on an assumed PNLF. IBuzzoni et al.l (2006) 
collated and presented a for both local group galaxies 
and a sample of early type galaxies (many of which we 
have plotted in figure [3]). This work, and the references 
therein, demonstrates that there are significant varia- 
tions in a. While some galaxies have a in the range of 
the upper limit, discussed in section |4~T1 others are a fac- 
tor ^7 lower. Despite these variations, most galaxies are 
found to have a in the range of 1 — 6 x 10~ 7 PN/L^ (e.g. 
i Jacobv fc Ciardulld[l99l IBuzzoni et alJl2006t [Ciardullol 
12010ft . 

Studies of the PN population in galaxies have ob- 
served correlations between a galaxy's properties and 
its population of PNe. Although the reason is yet to 
be fully understood, it is observed that more massive, 
metal rich galaxies tend to ha v e have lower a (e.g . 
Peimbertl [19901: IHui et all fl993t ICiardullo et all 120051 : 
Buzzoni et al.ll2006l ). Another correlation is observed be- 
tween a and the far-ultraviolet (FUV) excess of a galaxy 
(its FUV-V color). FUV bright gal axies appear to pro- 
duce less PNe per solar mass (e.g. Uacoby fc Ciardullol 
[19931 ICiardullo et al.l 12001 IBuzzoni et al.ll2006D . Such 
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an effect may be expected. The FUV excess in early 
type galaxies is now thought to be due to extreme 
horizontal branch ( EHB, or subdw arf B) stars (e.g. 
iDorman et all 119951: lO'Connelll 11999( 1 . Such stars may 
skip the AGB phase of stellar evolution entirely and pro- 
ceed directly to the white dwa rf phase as AGB-manque 
stars (jGreggio fc RenzinHll99ClD . Therefore, stellar pop- 
ulations which produce relatively more EHB stars are 
likely to be under abundant in PNe. It should be noted 
that the masses, metallicities and ages of early type 
galaxies are related - with massive, metal rich galaxie s 
having older stellar populations (e.g. lTrager et al.1 2000). 
Therefore, any/ all of these paramters could drive the 
observed PN relationship. 

4.2. Known PNe in globular clusters 

In table [5J we list all of the GC PNe that are cur- 
rently proposed. In th e Galactic GC syst em only four 
PNe have been found (jJacobv et al.1 119971 ). For extra- 
galactic GC systems relatively few PNe are known, par- 
tially due to a lack of data and studies aimed at identi- 
fying emission lines in GCs. Three PNe have be en pro- 
posed in extragalact ic GCs in: the Fornax dSph. | Larsen 
2008); NGC 3379, IBergond et all (l200fih . IPierce et al " 



2006); and NGC 7457. IChomiuk et all (|20081. In addi- 



tion to these GC PNe, [OIII] emission has been observed 
from several other clusters. One cluster in M 31, ap- 
pears to have [OIII] emission, which may be due to a 
PN, but this is yet to b e studied (as can be seen in the 
spect ra of the GC B115: IKim et"ld1l2lM ICaldwell et alJ 
12011(1 . Emission observed from several other extragalac- 
tic GCs may not be be due to PNe. These include two 
clusters in M87 (where the OIII emission is thought to 
originate fro m background fila ments in the central re- 
gion of M87; iCohen et al1ll998t l. a cluster in NGC 5128 
(which was identified as the first PN in an extragalac- 
tic GC, but higher resolution spectroscopy has subse- 
quently shown that the OIII emission is double peaked, 
suggesting velocities which exclude a PN as the origin; 
iMinniti fc Rejkuball200l iPeng et al.ll2004( l and a second 
source in NGC 7457 (with an un certain origin, but po s- 
sibly from a supernova remnant: IChomiuk et aLl l2008). 

Table [5] lists the properties of these PNe and the clus- 
ters which host them. It can be seen that the PNe span a 
wide range of [OIII] luminosities, from close to the tip of 
the PNLF to 9 mag fainter. Because of the low number 
of GC PNe known, and the varying survey depths, the 
GC PNLF is hard to constrain from this sample. How- 
ever, there is no evidence for either brighter or fainter 
PNe being favored in GCs, relative to Galactic regions. 
It is also interesting to consider the properties of the clus- 
ters which host a PN, since this may provide clues to the 
origin of GC PNe. The left panel of figure [21 shows the 
metallicity and luminosity of the Galactic GCs (crosses) 
and of the GCs which host PNe (red squares). As dis- 
cussed in section PL 1.31 a correlation is observed between 
a and the metallicity of early type galaxies, with metal 
poor galaxies hosting more PNe. It can be seen that 
no metallicity effects are observed in these GCs. A K- 
S test confirms that the metallicity of the PN hosting 
GCs is consistent with being drawn randomly from the 
Galactic GC population. The PN hosting clusters do 
favor more massive clusters. This may simply be a con- 
sequence of the greater number of stars in these clusters. 



Uacobv et"aLl (|1997( ) also noted that 2/4 of the Galac- 
tic GCs, which host a PN, also host an LMXB. From 
this observation, they suggest that PNe and LMXBs 
may form preferentially in the same clusters (albeit with 
a significance < 3cr). Such a correlation would be in- 
teresting, since the presence of an LMXB in a cluster 
is considered indicative of the formation of tight binary 
systems. Therefore, this may suggest a possible binary 
origin for the GC PNe. A more direct way to study 
such a relationship is by considering the stellar interac- 
tion rate in the dense cores of these clusters, T c cx p 3 / 2 /r 2 
(whe re p c and r c are the c luster's core density and ra- 
dius; ^rbunT^Hullilp). If the formation of PNe is 
related to binary interactions, then they should be found 
primarily in clusters with high r c . A strong correlation 
is observed between T c and the f ormation of LMXBs i n 
the GC sys tems of the MW (e .g. IVerbunt fc rbul[T987l) 
NGC 5128 ([Jordan et al.ll2007l) and M31 (|Peacock et al.1 
I2010bf) . The right panel of figure [5] shows the stellar in- 
teraction rate for all Galactic GCs and those which host 
PNe (red squares) and LMXBs (blue circles). Both the 
LMXB and PN hosting clusters have higher than aver- 
age r c . However, due to the low number of PNe, the 
difference between the Galactic GCs and the PN GCs is 
not significant. Furthermore, it can be seen that there is 
a strong correlation between the luminosity of a cluster 
and its interaction rate. While the LMXB hosting clus- 
ters have significantly higher than average r c for their 
luminosity, strongly supporting a dynamical formation 
origin, the PN GCs do not. 

From the known GC PNe, we therefore conclude that 
there is no evidence for enhanced dynamical forma- 
tion of PNe in these GCs. The partial correlation be- 
tw een the presence of a PN and an LMXB, proposed 
by Uacobv et al.1 JT997), is likely a consequence of the 
high stellar mass in the two clusters which host both a 
PN and an LMXB. The PNe are found preferentially in 
more massive clusters. This is consistent with their for- 
mation from either isolated stars or primordial (rather 
than dynamically enhanced) binary interactions. 



4.3. a in globular clusters 

Due to the low number of GC PNe currently known, 
and the lack of systematic searches for such objects in 
large GC samples, a in GCs is currently poorly con- 
strained. An empirical estimate for a can be produced 
based on the Galactic GCs, al though this estimate suffers 
from small number statistics. Uacobv et al.1 (|1997t ) identi- 
fied four PNe in Galactic GCs with a combined stellar lu- 
minosity of 2.4xlO 7 L , implying a of 1.7x 10~ 7 PN/L Q . 
This survey covered the brightest ^9 mag of the PNLF, 
significantly deeper than most extragalactic surveys of 
PNe can achieve. For surveys sensitive to the brightest 
2.5 mags of the PNLF (such as the data considered in this 
paper), only one PN, that in Pal 6, would be detected in 
the Galactic GCs (see table [2]). It should also be noted 
that the association of this PN, with Pal 6, is relatively 
uncer tain (with a 0-5% chance that they are uncorre- 
lated. ljacoby et aT1ll997t ). The absence of any PNe in the 
NGC 4472 cluster data provides further important con- 
straints on the formation of PNe in GCs. The combined 
stellar luminosity of these clusters is ^2xl0 8 L Qi v (~8x 
that of the Galactic clusters) . This corresponds to a bolo- 
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TABLE 2 

Planetary nebulae in globular clusters 



Names 


(m-M) 1 


My 


[Fc/H] 1 


log(rc) 1 - 2 


F5007 


1115007 


L5007 


M5007 


M*-M 


references 


host galaxy/GC/PN 










erg/s/cm 2 




X 10 35 crg/s 








MW/M15/Psl-K648 


15.09 


-9.19 


-2.37 


6.83 


1.87E-12 


15.58 


0.242 


0.49 


4.97 


Adams et al. (1984) 


MW/M22/GJJC1 


12.53 


-8.50 


-1.70 


5.63 


3.0E-13 


17.57 


0.004 


5.04 


9.52 


Gillett et al. (19891 


MW/Pal6/JaFul 


13.82 


-6.79 


-0.91 


5.28 


1.25E-10 


11.02 


5.032 


-2.8 


1.68 


Jacobv et al. (1997) 


MW/NGC6441/JaFu2 


15.32 


-9.63 


-0.46 


7.17 


1.9E-13 


18.06 


0.031 


2.74 


7.22 


Jacobv et al. (1997) 


NGC3379/gc771/PN 


29.91 


-8.65 


-1.40 




1.54E-17 


28.29 


1.700 3 


-1.62 


2.86 


B06, P06, Pll 3 


Fornax-dSph /H5/PN 


20.68 


-7.40 


-1.73 




5.0E-15 


22.01 


0.112 


1.33 


5.81 


Larsen (2008) 


NGC7457/GC7/PN 


30.55 


-8.20 


-0.40 




6.04E-17 


26.81 


12.000 


-3.74 


0.74 


Chomiuk et al. (2008) 



1 for the Galactic GCs, this data was obtained from the Harris catalog (2010 edition; Harris 1996) 
2 the stellar interaction rate in the cluster core, r c oc p 3 / 2 /r 2 (see e.g. IVerbunt &: HutH198^) 

3 this source was identified by both Bcrgond et al. (2006, B06) and Pierce et al. (2006, P06). However, the [OIII] luminosity is not quoted 
in these papers. We measured this from the Flames/Giraffe spectrum used by Bcrgond et al. (2006). This spectrum was obtained and 
reduced in a similar fashion to the NGC 4472 data, discussed here, and the line measured using IRAF/SPLOT. 
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Fig. 2. — Absolute V-band magnitude of Milky Way GCs (crosses) as a function of [Fe/H] and stellar interaction rate in the cluster core, 
r c . The PNe hosting clusters are highlighted as red squares. For the magnitude and metallicity, we show all known GC PN, but r c is only 
available for the four Galactic PN GCs. For comparison, the LMXB hosting clusters are highlighted as blue circles. 



metric luminosity, Lboi ~2.8xl0 8 L Q (assuming L DO i for 
an ol d simple stellar population is ~1.4xLv. iMarastonl 
2005). The detection of no PNe in these GCs therefore 
corresponds to «2.5 <0.8x 10 _8 PN/L Q (where the upper 
limit is b ased on the 90% confidence level for a non- 
detection, iKraft et al.l ll991). This is smaller than that 
measured for the similar stellar populations of the Galac- 
tic GCs. However, the numbers are consistent, given the 
low number of sources detected. 

Both GC systems host significantly less PNe than po- 
tentially poss i ble fo r their stellar mass, see section I4.ll 
Uacobv et al.l (|1997l ) have previously demonstrated that 
the Galactic GCs host less PNe than this limit, with a 
significance of 3.1a. The non-detection of any PNe in the 
larger stellar mass covered by NGC 4472 's GCs strength- 
ens this finding, with a < 7x lower than predicted. This 
suggests that, either some of the AGB stars in GCs do 
not evolve into PNe, or the lifetimes of PNe in GCs are 
significantly shorter than the dynamical timescale. The 
lower a in these GCs could be due to their old stellar 
populations having PN central star masses that are too 
low to produce observable PNe. However, given the sim- 
ple stellar populations and similar (old) ages of the GCs, 



this raises the further question of how the PNe that are 
observed in the Galactic GCs were produced. One possi- 
ble explanation is that the GC PNe are the evolved stage 
of blue straggler stars. Blue stragglers are known to be 
present in GCs and can have masses significantly greater 
than the main sequence turnoff. If only blue stragglers 
are capable of producing PNe in GCs t hen their rarity 
may e xplain the lower rates observed. iCiardullo et al.l 
(2005) have previously discussed the possible role of blue 
stragglers in producing the brightest PNe observed in 
galaxies. Unfortunately, the low number of PNe in the 
Galactic GCs makes a direct comparison with their blue 
straggler population impossible. However, observations 
of the Galactic GC PNe lends some support to this the- 
ory. The masses of the central stars in two o f these PNe 
have been measured at 0.62 and .75M W ()Alves et al.l 
120001 Harrington fc Paltogloul 119931 ). Such stars would 
have main sequence masses much higher than the main 
sequence turnoff mass in these clusters. 

As discussed in section 14.1.21 binary interactions may 
be required to produce PNe, with only 13% of stars 
thought to have binary companions in the required range. 
Therefore, if binary interactions are required to pro- 
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Fig. 3. — Number of PNe per bolometric luminosity, a, as a function of the host stellar population's metallicity ([Z/H], left) and velocity 
dispersion (cr, right ). Crosses show a for a sample of early type galaxies, the MW bulge, LMC and SMC (taken from tables 4 and 6 of 
IBuzzon i ct al. 2006<). Filled squares show the MW's and NGC 4472's GC systems. See text for details on the data used. The dashed line 
is the predicted a assuming all AGB stars produce PNe with with lifetimes close to the dynamical timescale of the nebula. It can be seen 
that the observed a in these galaxies increases with decreasing mass and metallicity. However, this correlation does not appear to extend 
to the lower mass and metallicity GC systems. 

duce PNe, one would expect a corresponding reduction 
in a. The resulting rate would be in reasonable agree- 
ment with both the MW's GCs and the upper limit im- 
posed by NGC 4472's GCs. However, this does not ex- 
plain the lower fraction of PNe in GCs compared with 
some galactic populations. A potentially relevant ob- 
servation is that primordial binary fractions may be 
lower in the Galactic GCs than in the Galactic field (al- 
though such fract i ons a re very difficult to determine; e.g. 
IMoni Bidin"et~aTI 12001 ) . If binary interactions are re- 
quired, then differing binary fractions/ parameters may 
correspond to the different as observed. 

It should be noted that, binary systems are signifi- 
cantly effected by the dense stellar environments in GCs. 
New binary system s can be formed vi a tidal capture 
(iFabian et al.lfl975f ) or direct collisions (jVerbunt k, Hutl 
an d binaries are driven to shor ter periods via dy- 
namical (|Heggidll975| ) and exchange ([Hills fc Davlll976( ) 
interactions. The result is that GCs contain more close 
binaries, per stellar mass, than galactic field regions (as 
confirmed via the dr amatic exces s of accreting compact 
objects in GCs, e.g. iClarkl rT"975f ) . Therefore, if binary 
interactions are important in producing PNe, one may 
expect this process to be enhanced in GCs. This does 
not appear to be the case in these GCs, which have a 
lower than that of many galaxies. This suggests that, 
close interacting binaries may not enhance the formation 
of PNe. In fact, by driving binaries to shorter periods, 
it may be that binary interactions reduce the number 
of PNe in GCs. If the binaries have shorter periods, 
and hence their stars interact in the early AGB phase, 
or even during the RGB phase, the envelopes of these 
stars may be ejected from the system before the central 
star is hot enough to ionize the gas and hence produce 
a PN. Interactions with a companion star on the RGB 
have previously been proposed t o enhance the forma- 
tion of EHB stars in dense GCs (jBuonanno eTaTl [19971 : 



lHan et al.H"2007HPeacock et al.ll2010aD . Since EHB stars 
may skip the AGB phase entirely, such interactions may 
reduce the number of PNe in these clusters. A similar 
argument may be related to the low a of FUV bright 
ellipticals, which are thought to have higher fractions of 
EHB stars. 

Figure [3] compares a for the MW's GC system, 
NGC 4472's GC system and a sample of gal axies (taken 
from tables 4 and 6 of IBuzzoni et al.l 120061) . It can be 
seen that a for the Galactic GCs is consistent with that 
measured in many early type galaxies. However, there 
is a large uncertainty in a for the Galactic GCs, due 
to their low combined stellar mass and the low number 
of PNe associated with them. The non-detection of any 
PNe in NGC 4472's clusters requires that these GCs have 
one of the lowest as of any stellar population observed 
(only 18% of the galaxies shown in figure [3] have a be- 
neath the upper limit required for these GCs). However, 
some galaxies do have a that is consistent with the limit 
suggested for NGC 4472's GCs. 

In figure [3] we plot a as a function of the metallic- 
ity and velocity dispersion (cr) of the stellar populations. 
For the early type galaxies their velocity dispersion and 
Mg2 l ick indices are taken from table 4 of IBuzzoni et al.l 
(2006). We convert Mg 2 to [Z/H] using stellar popula- 
tion model predicti ons for a 12Gyr old population, as- 
suming [q/H]=0.3 (j Th pmas et alj |200l iMarastonl 120051: 
iThomas et all UMljlZ]. For local group galaxies the 
m etallicities were taken f rom the RGB star obse r vations 
of ICarrera et~aTl (l2008aL SMG) ICarrera et al.l (I2008bl 
LMC) and iMcConnachie et al.l ([20051, NGC 205). The 
metallicity and velocity dispersion o f the G alact ic bulge 
were taken from iMcWilliam fc Rich! (|1994| ) and iMinnitil 
(1996), respectively. Similarly to individual stars in 
galaxies, the metallicity of individual GCs in both the 

Iv taken from http:/ /research. icg.port.ac.uk/~thomasd/tmj. html 
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MW's and NGC 4472's GC systems varies (with GC val- 
ues in the range -2.5<[Z/H]<0). We plot the average 
values of the GCs in both systems. The metallicity and 
a quoted for the MW GC system is th e mean of th e 
GC values, as presented in the catalog of IHarrisI (|l996). 
For NGC 4472's GC system we have converted the GCs 
mean B-R color and V-band m agnitudes (taken from the 
catalog of iRhode fc Zepj|200l to [Z/H] and tr using the 
correlations observed between these parameters in the 
MW's GCs. The previously observed correlation between 
a and the mass/ metallicity of the galaxies can be seen. 
However, these correlations do not seem to extend to the 
relatively low mass, metal poor, GCs - which have much 
lower a than would be predicted by the galaxy relation- 
ships. The reason for these correlations in the galaxy 
data is still uncertain, making it is hard to understand 
the implications of them breaking down in GCs. If the 
PN population is effected directly by either metallicity 
or mass, the lower rates observed in GCs could be due 
to other reductive processes, such as an age cut off or 
binary interactions, dominating over these effects. An- 
other potential explanation for the observed rates is that 
the correlations observed are driven by the ages of the 
stellar populations. The massive, metal rich, early type 
galaxies will t end to have older ages, closer to those of 
the GCs (e.g. iTrager eFaT][2000l) . Therefore, if the cor- 
relations observed in figure [3] are due to population age, 
rather than metallicity or mass directly, one may expect 
the GCs to have a similar PN population to the massive 
high metallicity older galaxies, as observed. 

Investigating age effects directly in this sample of GCs 
is hard due to their similar (old ] ages. Our results a gree 
with the low values pre dicted bvlBuzzoni et all (|2006() for 
these old populations. IBuzzoni et all ( 20061 ) also predict 
that young stellar populations should have low a. This is 
because the young central stars in the PNe evolve quickly, 
shorte ning the lifetimes of the PNe. lLarsen fe Richtlen 
(2006) have previously studied a sample of 80 massive 
young clusters (with ages around 100 Myrs) in several 
nearby galaxies. They identified four PNe candidates 
in these clusters. This suggests that these young popula- 
tions do have a low a of 0.6 x 10 _7 PN/L© - similar to the 
low value found for the MW's and NGC 4472's GCs (with 
ex pected ages >10 Gyr s). As can be seen from figure 15 
of IBuzzoni et all (|2006f) . higher values of a would be ex- 
pected at intermediate ages of a few Gyrs. Investigating 
the PN population of such clusters would be highly de- 
sirable. Unfortunately, identifying int ermediate age clus- 
ters has proved very challengi ng (e.g. lKundu et al.ll2005l : 
lLarsen et all 120051 IZeBI 120091) . 

5. BLACK HOLE LMXBS IN GLOBULAR CLUSTERS 

5.1. Nebular emission from ultraluminous LMXBs 

In section [3J we demonstrated that the only cluster in 
our sample which shows [OIII] emission is, the black hole 
(BH) LMXB hosting cluster, RZ 2109. This cluster is in- 
teresting because it hosts one of the few known BHs in a 
GC. One of the most conclusive signatures of a BH pri- 
mary in an LMXB, is to detect variable X-ray emission 
that greatly exceeds the Eddington limit for an accret- 
ing neutron star. Using this method, B Hs have been pro- 
posed in this GC ([Maccarone et all2007 ) , a second GC in 
NGC 4472 (not covered by our survey. iMaccarone et all 



l20Tfl) . a GC in NGC 3379 (IBrassington et alllMlOh and 
in tw o GCs in NGC 1399 (Shi h et all 120101: llrwin eTall 
2010). It is thought that LMXBs accreting at such high 
rates (at or above the Eddington limit) can drive strong 
winds from the syste ms, which will likely be ionized 
by th e central LMXB (jBegelman et all 120061: IZepf et all 
2008). This would suggest that ultraluminous LMXBs 
may generally be associated with nebula emission and 
is the proposed source of the [OIII] emission observed 
in RZ 2109. Unfortunately, the other BH-LMXB candi- 
dates cannot be used to test this ubiquity. The [OIII] 
emission from three of the candidates is currently un- 
known. The other ultraluminous BH candidate is as- 
sociated with [OIII] and [Nil] emission lines. However, 
this source has different properties to RZ 2109 and may 
be an intermediate mas s BH accreting f rom a tidally 
disrupted while dwarf ([Irwin or a stellar 

mass BH ionizing the wind of a R Corona Borealis star 
(M accaronc fc Warnenl20Tl . 

5.2. Ultraluminous LMXBs: BHs or beamed NSs? 

The detection of X-ray emission, with luminosities in 
excess of the Eddington limit for an accreting NS, pro- 
vides some of the best evidence for BHs in GCs. How- 
ever, it has also been proposed that such luminosities can 
be produced by an accreting NS in an ultracompact X- 
ray binary, i f the X-ray emissi on is geometrically beamed 
towards us (|King et al.l 120011) . The Eddington limit for 
accretion on to a 1.4M Q NS is ~ 1.8 x 10 38 erg/s. If 
we assume hydrogen poor accretion, as may be expected 
from a white dwarf donor in an ultracompact LMXB, 
slightly higher Eddington limits of ~ 3.5 x 10 38 erg/s can 
be reached. The observed luminosity (L D b s ) can then be 
enhanced, along the line of sight, via geometric beam- 
ing, to give L bs=LEdd/&- Here, b is the beaming fac- 
tor (—il/An). Hence beaming factors of 0.1-0.3 would 
be required to produce the observed luminosities of the 
BH LMXB candidates. Significantly, in such models, the 
[OIII] emission should not be beamed. This is because 
the nebula emission has to come from much larger scales 
than the X-ray emitting region of the binary. 

If [OIII] emission, such as that observed in RZ 2109, 
is associated with LMXBs accreting at super-Eddington 
rates, this provides a way to discriminate between 
beamed-NS and unbeamed-BH models. Under the beam- 
ing scenario, there should be, on average, 1/6 additional 
LMXBs accreting at super-Eddington rates, but beamed 
away from our line of sight. If these systems are present, 
their (unbeamed) [OIII] emission should be similar to 
that observed from binaries in which the X-ray emission 
is beamed towards us. We would therefore expect to ob- 
serve [OIII] emission from (1/6) x iV S Edd clusters, where 
^VsEdd is the number of LMXBs with observed super- 
Eddington accretion rates. From the spectra studied 
in this paper, we can confirm that there are no addi- 
tional RZ 2109 like sources in this sample. Based on 
the observation of one bright LMXB in these data, the 
above beaming factors would predict an additional 3-10 
[OIII] sources. Our non-detection of any sources there- 
fore favors mild beaming factors >0.4 (based on a con- 
fidence level of 0.90 for the non-detection). It should 
be noted that, this conclusion assumes that the one 
super-Eddington LMXB in our sample is representative 
of the global population. Currently, the number of super- 
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Eddington LMXBs is quite poorly constrained, since we 
are dealing with very small numbers of sources. Con- 
sidering this uncertainty, our data suggests that high 
beaming factors are unlike ly, bu t milde r factors of ^0.3, 
such as those proposed by iKingl (|2011[ ). cannot be ruled 
out. Therefore, we conclude that the non-detection 
of additional [OIII] sources in these clusters favors the 
BH LMXB hypothesis with mild or no beaming. We 
stress that more [OIII] data and new detections of super- 
Eddington LMXBs should help to better constrain these 
conclusions in the future. 

Additional to the lack of [OIII] emission in other 
GCs, there are difficulties in reaching the high X-ray 
luminosities using the beamed NS LMXB hypothesis. 
Kind (|201lD suggests that the BH candidate in RZ 2109 
(XMMU 122939.7+075333), another NGC 4472 GC BH 
candidate (CXOU 1229410+075744) and the hyperlumi- 
nous source HLX-1 can be explained as super-Eddington 
NSs in ultracompact LMXBs. The reason for invoking 
ultracompact LMXBs is the claim that they can have 
higher effective Eddington luminosities, so that they need 
to exceed their Eddington luminosities by smaller factors 
in order to appear at the observed luminosities. It is 
widely discussed that for Type I X-ray bursts, this is the 
case - as discussed above, hydrogen-poor material will 
have an effective Eddington limit almos t twice as high 
as th at for hydrogen rich material (e.g. iKuulkers et al.l 
2003). However, this is not the case for ultraluminous 
X-ray sources, which have characteristic temperatures 
nearly an order of magnitude lower than those for Type 
I X-ray bursts. At kT ~ 0.1 - 0.2 keV, as is typically 
seen for ULXs, the effective Eddington luminosity for 
material from a carbon-oxygen white dwarf becomes dra- 
matically lower than it is for hydrogen rich material, ex- 
cept at very low densities. This is because the opacity 
stops being dominated by Thomson scattering, and be- 
gins ^to_j2e_d^mhiated bycarbon and oxygen edges (see 
e.g. Ilglesias fc Rogers! 1 1993] ). 

If we take parameter values suggested by IKingl ([201 ID 
for XMMU 122939.7+075333 (in RZ 2109); an accretion 
rate of 10 19 g/sec and a 1.4M© NS. The spherization ra- 
dius, where the accretion disk solution becomes locally 
Eddington limited, occurs at 4 x 10 7 cm. The accretion 
disk solution should approximate a slim disk not too far 
outside this radius. One can estimate the density of the 
accretion flow as Mt V i sc /(4TrRg P ), where M is the accre- 
tion rate into the outer disk, t V i SC is the viscous timescale 
from the spherization radius, and R sp is the spherization 
radius. Setting t visc to a~ j L s (H/R)~ 2 t^, where H is the 
scale height of the disk, a V i S is the viscosity parameter 
in standard accretion disk models, R is the radius of the 
accretion disk, and t$ is the orbital period in the accre- 
tion disk and taking (H/R) = 0.1 and a V i S — 0.1 (both 
typical parameters), we find that the density of the ac- 
cretion disk will be about 10~ 3 g/cm 3 . On e can then ex- 
amine figure 9 of Ilglesias fc Rogers! ([1993), and find that 
the opacity per unit mass will be about 10 times as large 
for a carbon/oxygen mixture at 2-3 million K than it will 
be for pure hydrogen gas at the same temperature. As 
a result, the effective Eddington limit will be a factor of 
10 lower than the standard Eddington luminosity, rather 
than a factor of 2 higher than the standard Eddington lu- 
minosity, and both a larger spherization radius, by a fac- 



tor of about 3 and a larger effect from geometric beaming 
(by a factor of about 1 0) are needed relative to the pro- 
posal from (|Kind l20Tll ) . The accretion disk would then 
be expected to peak in the extreme ultraviolet, rather 
than in the X-rays, and the model becomes incapable of 
explaining the high X-ray luminosities observed. If we 
consider a BH, rather than a NS accretor, the density 
of the gas will be sufficiently low that it will be highly 
ionized - enough that the opacity will not be far different 
from the Thomson opacity. The result is that, for a BH 
accretor, there is only a minor perturbation to the ac- 
cretion solution by having a donor made of carbon and 
oxygen, rather than a donor of roughly solar composi- 
tion. 

Additionally , for the other BH candidates considered 
bv lKmel (f201lb . CXOU 1229410+075744 and HLX-1, the 
model fails to fit their observed X-ray spectra. In both 
cases, the observed thermal components are far too hot 
to be explained by a model involving super-Eddington 
accretion without the model having much larger beam- 
ing fa ctors than those invoked in the model of IKingl 
(poll . For CXOU 1229410+075744, the X-ray spectra 
are well-fitted by disk b lackbody models with kT ^1 keV 
(jMaccarone et alJl^Ml . The model of EEi (poTTI ) re- 
quires an accretion rate about 11 times the accretion rate 
required for the source to be at the Eddington luminos- 
ity, and hence predicts a characteristic temperature for 
the source's thermal c omponent of about 0.1-0.2 keV. For 
HLX-1, the model of En3 (poll requires an accretion 
rate 170 times that needed for the source to be at Ed- 
dington, meaning that the thermal component should be 
in the extreme ultraviolet, undetectable by X-ray obser- 
vatories, rather th an at about 0.2 keV as observed (e.g. 
iFarrell et al.ll2009]) . 

6. CONCLUSIONS 

Our survey of 174 massive GCs in NGC 4472 iden- 
tified only one cluster with [OIII] emission. As dis- 
cussed above, this source is thought to be associated 
with the LMXB in this cluster, which is accreting at 
super-Eddington rates. The high X-ray luminosity of 
this source, combined with the absence of any similar 
[OIII] sources in the other GCs, favors the theory that 
this emission is from a BH LMXB, rather than a NS 
LMXB that is geometrically beamed. The absence of 
any narrow [OIII] lines in these data also sets constraints 
on the formation of PNe in these clusters. This limit 
requires that GCs have some of the lowest PNe forma- 
tion rates observed. However, galaxies with similarly low 
rates are observed. We note that correlations, observed 
in elliptical galaxies between their PN population and 
the galaxies mass and metallicity, do not extend to the 
low mass and low metallicity region occupied by GCs. 
Potentially, this may point to other parameters, such as 
age, being the cause of these correlations. The formation 
rate of PNe in GCs is significantly lower than the num- 
ber of post-AGB stars predicted from stellar evolution, 
suggesting that some of these stars do not produce PNe. 
The lower rates observed are consistent with the lower 
rates predicted, if PNe require a binary origin. However, 
the low a in both the MW's and NGC 4472's cluster 
systems, suggests that the formation of PNe is not en- 
hanced in GCs due to dynamical formation of tight bina- 
ries. Indeed, dynamical interactions may act to reduce 
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the number of PNe, by enhancing mass loss on the RGB / 
early ABG branch. Larger samples of GCs with [OIII] 
data should be highly beneficial to better constrain these 
conclusions in the future. 
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